Quantum simulations based on measurements and feedback control 
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We propose a scheme for performing quantum simulations with atoms in cavities based on a photon 
detection feedback loop that requires only linear optical elements. Atoms can be stored individually 
without the need of directly interacting with one another. The scheme is able to simulate any 
time evolution that can be written as a sum of two-qubit Hamiltonians, .e.g., any next neighbor 
interaction on a lattice. It can also be made robust against photon losses. 

PACS numbers: 



The concept of building a quantum simulator goes back 
to the early days of quantum information theory, where 
Feynman 1] formulated the idea to build one special 
quantum device to simulate another one. This idea was 
forced by the impossibility to simulate any growing quan- 
tum system efficiently on a classical device because of the 
exponentially growth of the number of involved param- 
eters. However, this restriction can be ruled out by any 
quantum device due to the same exponentially behavior. 
The remaining problem is to find a well controlled quan- 
tum system which can be forced to evolve reliably under 
different time evolutions, especially under the physically 
mostly relevant next neighbor Hamiltonians. 

Several systems have been theoretically proposed as 
candidates for quantum simulation and enormous exper- 
imental effort has been put in the realization of the latter. 
Atoms are one of the most promising candidates for stor- 
ing a qubit, the basic quantum information unit for the 
quantum simulator. Trapped in optical lattices or ar- 
ranged in ions traps [1, Ql they allow for controllable next 
neighbor interactions that enable the simulation of e.g., 
the Hubbard-model Q or spin-lattices 0]. 

Quantum simulations are closely related to the gener- 
ation of entanglement since it is the presence of entangle- 
ment that prevents the system to be simulated by classi- 
cal algorithms. Generation of entanglement between e.g. 
two atoms is usually achieved by either using direct in- 
teractions between the atoms, or the interaction can be 
achieved indirectly by the exchange of photons . 

An alternative way to entangle two systems without 
ever getting into contact can be achieved by entangle- 
ment swapping. Two photons emitted by the two atoms 
can be jointly measured The outcome of the mea- 
surement defines a local feedback action one has to apply 
to finish the entangling protocol. This ansatz is avoiding 
the difficulties that arise when a single photon has to be 
re-absorbed by a single atom. Even in the case where 
we restrict to linear optical elements, entanglement can 
be reliable and efficient generated in a probabilistic way 
using measurement and local feedback action [lOj . 

Taking advantage of this idea for quantum simulation 
is hindered by two problems: (i) Using measurements to 
run a quantum simulation is an unusual concept since this 



(a) 



I'O) 





ll> 



FIG. 1: The atom stores a qubit in two ground states. In ad- 
dition several excited levels are used to create photons in the 
cavity. The double arrays indicate photon modes supported 
by the cavity, whereas the single arrows stand for external 
laser pulses to excite the atoms. In (a) only one photon mode 
is supported by the cavity. The photon qubit is encode is the 
presence or absence of a photon. In (b) two different photon 
modes are supported. 



does not lead to an unitary time evolution. The emitted 
photons have to be entangled with the atoms such that 
every measurement collapses the atom state to a random 
outcome. Thus, every single atom runs through its own 
random independent evolution and it may be hard to syn- 
chronize all the atoms to simulate one controlled unitary 
evolution, (ii) In a linear optical setup a complete Bell- 
measurement, as required for entanglement swapping, is 
not available In consequence, the desired interaction 
can only be generated in a non-deterministic way. In 
addition, a realistic photon is highly vulnerable to ab- 
sorption and measurement losses. 

In this letter we propose a scheme to circumvent those 
problems and to simulate effectively arbitrary next neigh- 
bor interactions among individually distance atoms with 
the help of emitted photons and measurements. Al- 
though each individual evolution is random, the sys- 
tem can be forced efficiently to follow any specific uni- 
tary next-neighbor evolution. The scheme can also be 
made resistant against any photon loss rates by keeping 
a backup copy of each photon. 

Photon generation: We consider atoms each in an opti- 
cal cavitywith the two ground states |0) and |1) realizing 
a qubit [7|. We assume that we can reliably apply any 
local unitary on the level structure of a single atom via 
lasers and external electromagnetic fields. To couple two 
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atoms we need the atoms to emit photons. For this pur- 
pose we need to include more internal levels of the atoms. 
We assume here two different possibilities depending on 
whether we want to encode the photon-qubit into two 
different polarizations or the photon occupation number 
(photon/ no photon) : 

(a) Photon/no photon: The atom has the two excited 
levels |eo) and |ei) , which decay via the same cavity mode 
into the ground state as shown in Figure [T] (a). The two 
ground states can simultaneously be changed in the two 
excited state by a laser pulse. By changing the strength of 
the laser pulse we can bring the ground states |0), |1) into 
any superposition of the ground and the excited states, 
i.e., |0) ^ VT^IO) + V^leo) and |1) ^ VT^|1) + 
y/e\ei). When the excited state decays the atom emits a 
photon such that we get the transformation 



|0) -^VT~s\0)\V)+V~e\l)\H), (1) 
|1) ^ VT~e\l)\V) + VemH) 



where \H) denotes a photon in the cavity mode, while | V^) 
stands for the vacuum, i.e., the absence of any photon. 

(II) Polarized photons: A similar kind of transforma- 
tion can be archived using the level structure showed in 
Figure [Ijb). The atom has four excited levels and the 
cavity supports two different cavity modes, called H and 
V. If the atom is in level |eo)(|ei)) it decays into |0) 
(|1)) and emits a ^-photon to the cavity, while when 
in state |eg) {\e[)) it decays to |1) (|0)) and emits a H- 
photon. With two laser pulses we can excite the ground 
states into a superposition of the four excited states 
|0) ^ VT^\eo)+V^\e'Q) and |1) ^ VT^\ei) + V^\e[), 
where the e can be controlled by the strength of the 
lasers. After the atom has decayed back into the ground 
states we get a transformation similar to ^ with the 
only difference that | V) stands now for a photon and not 
for the vacuum. 

In both cases the original qubit encoded in the ground 
state gets flipped with probability e where the flip can 
be recognized by a measurement of the emitted photon. 
This basic operation will be denoted as [4. 

Feedback loop: To get an effective interaction between 
two atoms we chose the following strategy: we apply Ug 
to the first atom and followed by a flip F to the 

second atom. The flip operations exchanges the states |0) 
and |1). The two (possible) photons emitted by the two 
atoms are sent onto a beam-splitter and measured. At 
one branch we add an extra phase of i to the iJ-photon 
(see Fig. EJ. 

If IV') denotes the initial state of the two atoms we will 




FIG. 2: Two distance atoms are effectively couplet by sending 
photons coming out of the cavity to a beam splitter. Depen- 
dent on the measurement outcome we apply a feedback loop 
on the atoms until the desired evolution has been done. 

get after emitting the photons the state 

C/e®FC/(i_,)|V') = (2) 

il~s)mVH) 
+ie(T (T\ip)\HV) 
+ ^/e{l-e)l(g)cr\i;)\VV) 
+i^ye{l^ e)a(g)l\tl;)\HH), 

where a is the first pauli matrix that flips the two 
states |0), After passing the beam-splitter the pho- 
tons are measured in the H, V basis. This correspond 
to an incomplete Bell-measurement of the two-photon 
pair, i.e., they are projected in one of the four states 
h) := \HV-VH),\+) := \HV + VH),\HH),\VV). 
The remaining two atoms are projected onto 

(i) K) (1— e)|V')— *£0'(8)cr|?/') with probability i((l — 

(u) |-|-) I— > {l~e)\ip)+iea^a\^) with probabihty ^((1 — 

(iu) \VV) 1-^ {1 (g) a)\ip) with probability e(l - e); 

(iv) \HH) 1^ (a(g) with probability £(1 - e). 

Let us first assume that the only possible outcomes 
would be |-|-), |— ). In both cases the resulting state can 
be written as e**'^*"^|V') = cos(t)|i/') + ism{t)a (g) 
with t equal to t = ± arcsin(e) « e. So we can simulate 
any time evolution with respect to the Hamiltonian a IS) a 
over a time t = ±e, but with a random time direction, 
i.e., by repeating we can simulate a random walk on the 
time axis of the evolution. If we aim for a time simula- 
tion V{t) = e**'^'^'^ we can repeat the protocol until the 
random walk by chance matches the desired time t. We 
can speedup this strategy by a feedback loop controlling 
the e parameter. In a first step we choose e in such a 
way, that sini = e giving us a chance of 1/2 to succeed 
with V{t) in the first step. We will assume t to be small 
such that we can approximate e = t va. the following dis- 
cussion, but in principle everything can be done exactly. 
In a case of a failure |— ) we end up with V{—t). By in- 
creasing £ to 2e we apply randomly V{±2t) in the next 
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step such that in the successful |+)-case we compensate 
the previous V{—t) rotation and in addition are left with 
the desired V{t) rotation. In case of further failures we 
each time double e and retry. This way we guarantee for 
a constant success probability of 1/2 in each round and 
we can simulate every evolution V{t) with exponential 
increasing probability in the number of rounds. 

Problems may occur if we measure \HH), \ VV) in be- 
tween. But in these cases we applied either l®a or a®! 
which can be corrected by local gates. Since our Hamil- 
tonian a® (J commutes with these unitaries we can ignore 
these measurement results and just correct the overall er- 
ror in the end. Note that this overall error only can be 
one of the three local unitaries 1® a,a ®l,a ® a. Since 
the probability of |— ) is at least 1/4 in each round 
the success probability is still exponentially good. So we 
can efficiently simulate V{t) for arbitrarily t up to some 
random but known Pauli-errors that commute with the 
evolution. 

Simulations for H — ® ui: Given V{t) and the 
possibility to do local unitary gates at each single atom, 
we can apply 



(1 + itcFk ® m) 



(3) 



for arbitrarily Pauli matrices ak^cru since V {t) = Uk'® 

uiV{t)u\ ® u\ with the local unitaries defined by Uk = 

1 _ 1 

(J^Gi ^. If the measurement outcome is \HH), \ VV) we 
pick up an error of 1 iK" cr; or cr^; (g) 1 which commutes with 
®. 

Simulation of an arbitrarily two-qubit Hamiltonian: 
Every two-qubit Hamiltonian can be written as H = 
J2ki ^kio-k fS) cr/. By successively applying V'^^eXa) we 
can simulate 

YlV^'iaXki) = l + ieY,^M<Jk(^<Ji + 0{e^) 

kl kl 

which equals e"^^^ up to first order in e. Note that each 
V'^''{eXki) can produce an error that does not commute 
with the following V'^ ' {sXk'i')- But since only products 
of Pauli-matrices are involved, the possible errors cither 
commute or anti-commute with the following evolutions. 
In the case that V'^ ' anti-commute only the time di- 
rection is inverted which we easily can compensate by 
exchanging the |+) and the |— ) case in above protocol. 
So we can simulate any product Yiki to a random 

but known error (Jr ® (Jp, i.e., \ip) ct,- (g) CTpJIfc; V^'^'IV')- 
This allows us to take advantage of the Trotter-formula 
111 to approximate e'*^ by 
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FIG. 3: To simulate a next neighbor Hamiltonian we can do 
m/2 of the rotations in parallel. We need to basic setup such 
that each cavity sends a photon on a beam splitter together 
with either its right or its left neighbor. 



Simulation of an arbitrarily next neighbor Hamilto- 
nian: In the same manner we can simulate any next- 
neighbor Hamiltonian on a lattice, e.g. the ID Hamilto- 



kl,x = l ^kl 



(x) (x) 
"k 



Ti^~^^\ where x labels the 



different qubits. Using the Trotter-formula and the V 
we can approximate e'*^ by 



kl 
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with a precession of the order of ^ up to an overall er- 
ror that is just a product of random but known Pauli- 
matrices. The same approximation can be made for any 
lattices with next-neighbor interactions or for any arbi- 
trary Hamiltonian consisting of sums of two-qubit Hamil- 
tonians. 

Time estimation: Assume we want to simulate a 
Hamiltonian H — X)"=i-^i' where denotes a two 
qubit Hamiltonian. Using the Trotter formula we can 
approximate the time evolution e*^* by applying Amn 
I^'^'-rotations. To get a constant precision we have to 
scale n at least linearly with m such that we need O(m^) 
applications of a single I^'^'-rotation for a proper sim- 
ulation. For every I/'^'-rotation we need one successful 
application of the feedback loop which succeeds with 
a probability of about ^ (or at least j ). If we ap- 
ply 2m^ -I- ^°^^* ?7i, steps we get at least successfully 
feedback loops with a probability that exponentially ap- 
proaches 100%, e.g. 97% for const — 3. Therefore we 
can simulate any Hamiltonian with m two-qubit terms 
in a time that scales quadratically with m. 

In a ID m-particle next neighbor setups the Hamil- 
tonian will consist out of m terms. Since only two 
qubits are involved in each local term we can speed up 
the simulation by applying at a time m/2 of the steps 
in parallel (see Fig. [3]). Assume we want to make 
m/2 gates in parallel in a time that allows for clog(m) 
feedback loops. Then each single gate succeeds with 
(1 - 2-^i°8('")) = (1 - m-'=i°s(2))^ gQ ^jjg ^Q^g^i proba- 
bility of success of the m/2 gates is (1 — m~'^'°s^^^)™/^ 
which approaches 1 as long as clog(2) > 1. This whole 
procedure has to be repeated 2m times to do the required 
m^ y'^'-rotations. Thus, we can simulate with arbitrar- 
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ily high probabihty m particles in a ID setup in a time 
scahng hke 0(rn log (m)). 

Errors: The protocol proposed so far is highly vul- 
nerable to photon losses and detection losses. A photon 
loss is especially fatal if we miss to identify a \HH) or 
\VV) measurement result, because we miss to identify a 
Pauli-error. 

In the case we work with polarized H and V photons 
we know whether we missed to detect a photon or not, 
since there are always two photons generated. We can 
use this fact to fight photon losses by the following idea. 
We will keep a backup-copy of the photon that we can 
use in the case that the real photon is lost. To store 
this copy assume a second 'atom' is present in each cav- 
ity. So we have two atoms (called A and B atom) with 
two internal levels each and we assume that we can ap- 
ply arbitrarily two qubit gates between them. In the 
first step the B-atoms will play now the role of the pho- 
tons to get the backup-copy before the photon is even 
generated: We start with the two atom state \iI>)aa' 
and the two extra atoms in state |00)bb/ and apply 
Ue now between the two atoms A and B (and A' and 
B'), i.e., \0)a\0)b 
|1)a|0)b 



(a) 



(b) 



VT^\0)a\0)b + V^\1)a\1)b and 



VT^|1)a|0)s + V£\0)a\1)b. After apply- 
ing this we are left with a state of form ([2]) where the 
photon modes V, H are replaced by the internal states 
|0), |1) of the two S-atoms. The next step is to copy the 
internal states of the B-atoms onto photons, such that 
we exactly get ^ and can continue the protocol. 

To this end we have to force the two i?-atoms to pro- 
duce a photon dependent on their internal level, i.e., 
|0)i3 |0)b|1^),|1)b Note that the B-atom 

will keep their quantum information and can be reused in 
the case of a photon loss. If we measure both _B-atoms in 
the |0 -I- 1), |0 — 1) basis then we destroy this information 
and are left with a state similar to ^ up to some random 
but known sign-phases that will change only the time di- 
rections. So we can use the above protocols up to the 
fact, that we have to include the measurement results of 
the i?-atoms to identify the direction of time evolution. 

Instead of measure the 'B' atoms in the |0 -f 1), |0 — 1) 
basis we can measure them individually in the standard 
basis |0), In this case \tp) stays unchanged up to one 
of the random but known errors 1 (8) cr;, (7^ iK" 1, cr^, (g) cr;. 
We do not get any time evolution, but are still able to 
continue the protocol. 

The main idea is now to first measure the photons and 
then the atoms. If both photons arrived we continue by 
measuring the _B-atoms in the |0-|- 1), |0— 1) basis. If one 
or both photons are missing we measure in the |0),|1) 
basis to protect our system against random unknown er- 
rors. Note that this measurement can be repeated until 
we successfully detected a photon. Doing so, we can han- 
dle any photon loss rates. 

Instead of putting two atoms in a cavity we can as well 
use one single atom with four internal levels. In the four 
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FIG. 4: (a) to fight photon losses the atom needs to store a 
second qubit to carry a security copy of the photon, (b) in 
the first attempt we try to establish an evolution by sending 
the photons on a beam-splitter. In case of success we erase 
the security copy, in case of failure we can reuse the copy by 
measuring in different basis. 



internal levels we can encode the A and the B atom. To 
this end we need an internal level structure that allows 
for storing two-qubits and the possibility to measure one 
of this qubits independently from the other. On possible 
level structure is shown in Fig. [4] The cavity supports 
two photon modes. If we excite |10), |00) to Icq), |eo) and 
|01) to \e[), |ei), we can measure the second qubit 
while the first qubit stays untouched. 

Quantum computation: Note that the same kind of 
setup can be used to do a quantum computation. If we 
apply V{t) to two atoms with t equal to tt/2 we get a 
two-qubit gate that is locally equivalent to a c-not gate. 
Together with local gates on single atoms this allows for 
universal quantum computation. This scheme is similar 
to the so called 'repeat until success'-schemes proposal in 
[l^ . The advantage we have here is that in our proposal a 
single two-qubit gate can be made photon loss resistant, 
whereas in 12] photon losses can only be attacked via 
the creating of cluster states which then are used as a 
resource for one way quantum computation. 

Conclusion: We have proposed a scheme for perform- 
ing quantum simulations based on a photon detection 
feedback loop that requires only linear optical elements. 
Atoms can be stored individually without the need to 
directly realize an interaction between two atoms. The 
scheme can simulate any time evolution that can be writ- 
ten as a sum of two-qubit Hamiltonian in a time scaling 
polynomially with the number of particles. Our scheme 
can be modified to be resistant to any amount of pho- 
ton losses by storing two qubit per cavity, whereas the 
second qubit carries a backup-copy of the photons. The 
same kind of setup can be used to do c-not gates and im- 
plement a universal quantum computer. So far the sim- 
ulation schemes needs several discrete round. One chal- 
lenge would be to make the interaction and the feedback 
control continuous and establish a continuous feed- 
back theory for a quantum simulator. We acknowledge 
support from DFG-Forschungsgruppe 635 and SCALA. 
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